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I.  INTRODUCTION 

Ion  beams  have  attracted  considerable  interest  in  recent  years  as  a  result 
of  their  potential  application  as  drivers  for  inertial  confinement  fusion. 

Due  to  the  state-of-the-art  of  accelerator  technology,  however,  theoretical 
advances  in  beam  propagation  have  focused  on  the  properties  of  high-energy 
(high-gamma)  electron  beams.  Some  important  simplifying  assumptions  are 
made  for  highly  relativistic  beams  which  do  not  apply  at  low-gamma.  In  past 
studies  (Refs.  1  and  2),  we  begarTto^^a^c  these  model  assumptions  to  address 
the  axi-symmetric  evolution  of  ion  beams.  In  this  SAI- -f-wefo  report  on  the 
propagation  of  low-gamma  ion  beams,  we  explore  the  resistive  hose  instability 
with  longitudinal  velocity-spread. 

For  energies  of  interest,  the  theoretical  models  of  electron  beam  propagation 
cannot  be  directly  applied  to  ions  because  of  their  low  relativistic  gamma  i 

(y)  factors.  For  example,  for  kinetic  energy  equal  to  50  MeV,y  *  99.  for  ^ 
electrons  compared  to  y  -  1.05  for  protons.  If  the  beam  propagates  in  the 
longitudinal  direction  (along  the  2-axis),  the  effective  inertial  mass  ,  l<4 
*  of  the  beam  particles  in  the  z-direction  is  given  approximately  bysy$,  where 
m  is  the  particle  rest  mass.  Thij  suggests  that  light  ions  are  considerably 
more  susceptible  to  longitudinal  acceleration  than  electrons  of  equivalent 
kinetic  energy;  the  longitudinal  particle-motion  cannot  be  ignored  as  in 
the  electron  case.  This  adds  a  new  dimension  to  ion  beam  propagation. 


'The  consequences  of  the  longitudinal  motion  on  beam  stability  can  be  divided 
into  three  categories: 


VT>  (1)  Longitudinal  spreading  of  the  beam  pulse  due  to  longitudinal 
velocity-spread  (longitudinal  temperature); 

(2)  Longitudinal  mass  (bunching)  instability;  and 

(3)  Resistive  hose  instability  with  longitudinal  velocity- 
spread. 
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Categories  (1)  and  (2)  are  unique  to  light-ion  propagation,  while  (3)  requires 
modification  to  existing  electron-beam  hose  stability  theory. 

In  order  to  address  issue  (1)  above,  we  developed  the  low-y  propagation  code 
LOGAP  (Ref.  1)  which  models  the  axi-symmetric  (monopole)  evolution  of  a 
single  pulse  (with  longitudinal  and  transverse  degrees  of  freedom).  LOGAP 
consists  of  four  basic  components:  electromagnetic  field  algorithm  GEM,  which 
generates  both  longitudinal  and  radial  field  profiles  in  real-time;  the 
highly-developed,  chemistry  package  BMCOND,  which  generates  the  conductivity 
profiles  in  the  ambient  gas  resulting  from  the  passage  of  the  beam  (Ref.  3); 
the  HIGAP  algorithm  for  solving  the  radial -envelope  equation  of  motion 
(Ref.  3  and  4);  and  the  longitudinal  dispersion  algorithm  LNGDSP  which 
describes  the  longitudinal  flow  of  beam  particles  via  a  fluid  model.  A 
detailed  discussion  of  LOGAP  is  found  in  References  1  and  2;  for  completeness, 
however,  we  provide  in  Appendix  A  excerpts  from  Reference  1  describing  both 
GEM  and  LNGDSP.  A  stability  analysis  of  the  GEM  field  algorithm  is  given  in 
Reference  2  and  is  reproduced  here  in  Appendix  B. 

Longitudinal  spreading  of  the  beam  pulse  is  a  concern  since  it  could  transform 
a  given  pulse  into  a  longer  one  with  lower  current  which  would  then  be  less 
hose-stable.  It  has  been  suggested  (Ref.  5)  that  the  spreading  might  be 
contained  by  the  self-induced  longitudinal  electric  field  at  the  front  and 
rear  of  the  pulse,  keeping  the  beam  particles  trapped  inside.  Preliminary 
results  from  LOGAP  to  date  have  failed  to  verify  the  existence  of  this  self- 
trapped  mode. 

It  should  be  emphasized  that  in  LOGAP  longitudinal  and  transverse  degrees 
of  freedom  are  essentially  decoupled.  If  we  partition  the  beam  particles  into 
subgroups  according  to  their  v  - velocity,  then  LOGAP  assumes  that  the 
current-density  of  all  subsets,  regardless  of  vz,  have  the  same  radial 
profile.  A  more  realistic  approach  would  be  to  allow  the  different  subsets 
to  expand  radially  at  different  rates;  the  low-velocity  subgroups  would 
expand  faster  since  their  particle  radial-motion  is  confined  by  a  relatively 


weaker  magnetic  pinch  force.  The  very  low  velocity  particles  (which  in  LOGAP 
contribute  to  the  longitudinal  spreading  of  the  pulse  tail)  may  "evaporate" 
in  the  radial  direction.  Therefore,  LOGAP  should  be  appropriately  modified 
to  account  for  transverse  evaporation  for  a  definitive  test  of  the  self- 
trapped  pulse  mode. 

The  longitudinal  bunching  instability  has  been  examined  by  Sloan  et.  at. 

(Ref.  6)  and  is  found  to  be  Landau  damped  when  the  longitudinal  rms  velocity 
spread  is  greater  than  some  minimum  value  [see  Eq.  (104)  below].  For  a 
^50  MeV,  10  ka)  proton  beam,  this  minimum  vales  is  %.075  v0>  where  \Q  is 
the  mean  velocity  of  beam  particles.  This  should  be  well  satisfied  at  the 
exit-port  of  an  Auto-Resonant  Accelerator  (ARA),  for  example  [see  Eq.  (101) 
below];  if  not,  the  bunching  instability  itself  and/or  the  self-induced 
electric  field  will  probably  ensure  that  the  rms  velocity-spread  within  the 
pulse  is  more  than  enough  to  stabilize  it. 

In  the  remainder  of  this  rep; .  we  examine  the  effects  of  longitudinal 
velocity-spread  on  the  resistive  hose  instability.  The  rigid  beam  model 
and  the  multi-disk,  distributed-mass  model  of  Lee  (Ref.  7)  are  both 
modified  to  account  for  the  phase-mixing  associated  with  the 
distribution  of  longitudinal  velocities.  Beam  particles  are  partitioned 
into  subgroups  according  to  their  vz~velocity  with  each  subset  experiencing 
an  independent  transverse  displacement.  The  resulting  equations  of  motion 
are  solved  for  mode  growth  in  the  context  of  an  initial-value  problem  with 
finite  pulse  length.  Upper-bounds  in  hose-mode  growth  are  found  via  saddle- 
point  analysis  as  well  as  conditions  for  "absolute"  hose  stability. 

The  results  of  this  preliminary  study  suggest  that  the  light-ion,  low-y 
beam  may  have  a  hose  stability  advantage  over  a  comparable  electron  beam. 

For  example,  the  (50  MeV,  10  ka)  proton  beam  with  longitudinal  rms  velocity- 
spread  of  %. 1  vQ  has  an  estimated  hose-stable  pulse  length  that  is  nearly 
twice  the  zero-spread  model  prediction.  With  an  rms  spread  of  .26  vQ,  it 
approaches  an  "absolute"  hose-stable  regime;  but,  if  the  pulse  spreads 
longitudinally  (i.e.,  no  self-trapping),  the  beam  current  decreases  until 
hose-growth  recovers. 


The  analysis  in  this  report  is  semi-quantitati ve  in  that  it  is  based  on  a 
single-mode  approach  with  various  simplifying  assumptions:  the  conductivity 
channel  is  assumed  fixed  in  space  independent  of  beam  evolution;  the  plasma 
return  currents  have  a  prescribed  radial-profile;  and  the  longitudinal  and 
transverse  equations  of  motion  are  decoupled.  In  a  comprehensive  treatment, 
these  assumptions  should  be  relaxed  and  the  equations  of  motion  solved 
numerical ly. 

In  Section  II,  the  basic  particle  and  electromagnetic  field  equations  are 
presented.  The  equilibrium  state  is  described  in  Section  III  with  emphasis 
on  the  isothermal  Bennett  distribution  and  the  longitudinal  velocity  distri¬ 
bution.  Hose-like  perturbations  are  considered  in  Section  IV;  the  rigid- 
beam  dispersion  equation  with  longitudinal  velocity-spread  is  derived  in 
Section  IV. A.;  and  the  modified  distributed-mass  model  and  the  single-mode 
dispersion  equation  is  derived  in  Section  IV. B.  The  dispersion  equations 
are  analyzed  in  Sections  V. A.  and  V.B.,  respectively.  And  in  Section  VI, 
we  interpret  the  results  in  the  context  of  the  ARA  and  longitudinal  mass 
instability. 


II. 


PARTICLE  AND  FIELD  EQUATIONS 


We  consider  a  beam  of  ions  with  mass  m  and  charge  propagating  with  a 
mean  velocity  vQ  in  the  positive  z-direction  through  an  ionized  gas.  The 
ambient  gas-plasma  is  characterized  in  the  present  study  by  a  real,  scalar 
conductivity  o  which  is  assumed  to  be  known  and  independent  of  beam  pertur¬ 
bations.  In  practice,  the  o-channel  could  be  strongly  affected  by  the 
passage  of  the  beam  through  the  gas  via  direct  collisional  ionization  and 
electromagnetic  field  breakdown,  but  these  effects  are  neglected  here. 

The  beam  current  1^  is  assumed  to  be  small  compared  to  the  Alfven  current 

dA): 


IL  « 


Y0$0  to10  2  - 


eZ.R, 
b  o 


=  Vo(t>  Tb  117  ka> 


(1) 


2-1/2 

where  yQ  =  (1  -  $Q)  '  is  the  relativistic  factor  for  the  mean  velocity 

vQ  =  Bqc,  me  is  the  electron  ..._ss  and  RQ  is  the  resistance  of  free  space 
(Rq  =  376.73  ohm).  Thus,  the  paraxial  approximation  is  adopted,  i.e., 


Ipll  «  po  ■  V*0  (2) 

where  |p^|  refers  to  the  amplitude  of  the  single-particle  transverse 
momentum. 


We  also  assume  that  a  is  large  enough  that  charge-neutralization  time  is 
small  compared  to  the  magnetic  decay  time,  and  hence,  the  displacement 
current  can  be  neglected.  Together  with  the  paraxial  assumption  (2), 
this  allows  a  reduction  of  Ampere's  law  to  the  form 


0  R  o  3A, 

VlAz  "  ~  aF  =  "°bz 


(3) 


where  the  electromagnetic  fields  are  adequately  described  by  the  longitudinal 
component  of  the  magnetic  potential  Az: 
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(4) 


ai  -  h  *  <az  ez> 

R  8A 

<5> 

The  transverse  components  (x,y)  are  denoted  by  the  symbol  (i),  and  Jfaz  is 
the  beam  current-density  in  the  z-direction. 

A  particular  beam  particle  with  velocity  v  =  (v^,vz)  satisfies  the  equations 
of  motion 


dp2  eZbR  f  3A 


.  r  r  3a 

r1 1  if  + 


dp,  eZ,  R 

JX  =  b  0  w  n  a 

dt  c  vz  V1  Az 


By  the  paraxial  assumption  (2),  we  ignore  the  variations  in  pz  and  assume 
that  v2  and  the  relativistic  factors  y  =  [1  -  (v2/c2)-  (v2/c2)]1//2  are 
constant  for  each  particle  orbit,  i.e., 


z  _  d 


=  dt  (Ymvz}  =  0 


The  mean  particle  momentum  is  expressed  by 


<p  >  ;  p„  =  y  jnv 
Z  O  0  0 


where  the  bracket  <  >  denotes  the  instantaneous  average  over  all  beam 
particles  in  a  disk  at  (z,t).  However,  it  is  also  assumed  that  pQ  be 
independent  of  (z,t)  over  the  time  scales  of  interest.  The  energy  of  a 
particle  with  vz  =  vQ  is  given  by 


E  -  y  me 
o  'o 


Opt*"'1 


It  is  convenient  to  express  the  particle  energy  and  momentum  in  terms  of 
the  momentum  variable  q  and  velocity  spread  v  where 


q  =  p  -  p  =  ymv  -  y  mv 
z  o  z  'o  o 


v  =  vz  -  v0 


(11) 


(12) 


with  q  «  p  .  By  (8),  both  q  and  v  are  considered  constants  of  motion.  To 
2  2U 

order  q  /(y^m^c^) ,  we  have 


2  2  Po^ 

E  =  ymc  -  y  me  +  +  -V  .  . 

!o  y  m  <v.3„  2yjn 


Pi 


'o'"  2y  m 
o 


1  0 


Pi  “  VV1 


(13) 

(14) 


Expanding  (11)  about  v  ,  it  is  easily  shown  that 


q  *  Yomv 


It  follows  that  approximation  (8)  is  valid  so  long  as 


lvll  vo 

<<  -x. 

y3  Y0 


(15) 


(16) 


Hence,  the  transverse  single  particle  equation  of  motion  (7)  can  be 
approximated  to  lowest  order  by 


dcr.  eZ.R 
1  bo 


v  V,  A, 


- —  _  -  y  V  |  r\ 

dt2  Tomc  0  1  z 


(17) 


Note  that  the  r.h.s.  of  Eq.  (17)  is  independent  of  v.  This  allows  us  to 
treat  particles  with  different  v  by  the  same  pinch-potential.  In  sub¬ 
sequent  sections,  beam  particles  are  divided  into  subgroups  according  to 
their  longitudinal  vz  velocity.  Eq.  (17)  implies  that  all  subsets,  regard¬ 
less  of  vz,  have  the  same  equilibrium  radial -profile. 
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*5  “3 

For  closure,  we  need  the  distribution  function  f(r,z,  v, ,  v  ,t)d"V  d^v  which 

JL  Z  r>  q 

represents  the  number  of  beam  ions  in  the  phase-space  volume  d^  dz  d  Vj_  dv2- 
The  beam  current  density  is  then  defined  by 

Jbz  =  eZb  /A.  A  vz  f  "  eZbV  (18) 


where  n  is  the  beam-particle  density.  From  continuity  and  Eqs.  (7)  and  (8), 
f  satisfies  the  Vlasov  equation 


if-  + 
at 


'•s.r  eZ,  R  v 

W  +  v  i},A  4  f  =  0 

1  -L  Z°Z  C  i  Z  PJ_ 


(19) 


It  is  useful  to  introduce  the  variable  £  defined  by 

S  5  vQt  -  z  (20) 

and  replace  the  independent  variables  (z,t)  by  (r,,z)  in  all  the  relevant 
equations.  The  derivatives  transform  as  follows: 


.  3_  9_ 

az  "  az  ac 


(21a ) 


(21b) 


d_  .. 
dt  _ 


z  az 


3  >V  f  E 

o  dz 


a_ 

az 


(21c) 


Employing  (20)  and  (21)  in  Eqs.  (3),  (17),  and  (19),  we  arrive  at  our  final 
set  of  model  equations  (dropping  subscript  z  from  now  on) 


*2  .  Roovo  aA 

* 1  A  "  W  =  *Jb 


dz 


2  rl 


eZ.R 
b  o 

Y  mcv  1 

'o  o 


Aa 


(22) 

(23) 


8 


f 


0,  -  cZ 

0 


r 


b  J  o  v±  jdv  V  f 

t 

b 

C — VlH'VPi' 


Af  eZ.  R 

dr  ,  ,  ,  b  o  o+  ,  i  *  ~  n 

vo  d?  +  V*lf  +  - r~ V,A^_  f  -  0 


(24) 

(25) 


These  are  employed  in  the  follov/ing 
instability.  We  begin  by  exploring 


sections  to  treat  the  resistive  hose 
the  axi- symmetric  equilibrium  state. 
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III. 


EQUILIBRIUM 


Our  equilibrium  beam  is  assumed  to  be  cylindrically  symmetric  about  the 
z-axis.  Adopting  cylindrical  coordinates  (r,0,z),  the  single  particle 
equation  (23)  becomes 


where 
which  by 


d2r  ?  h- 

dz2  * 


kB 


Wo 

Y  mv2 
o  o 


dr 


(26) 

(27) 


is  the  betatron  frequency  and  Aq  refers  to  the  equilibrium  potential 
(22)  and  (24)  satisfies 


IL  r  _ o  =  _  1  +  _ o  _ o  „  \  , 

r  dr  r  dr  Jbo  c  “m'  Jbo 

Jbo  =  eZb  1  v  fo 


(28) 

(29) 


In  (28)  we  have  assumed  that  the  plasma  return  current-density  is  represented 
by  -((m  Jbo  ,  where  the  current  neutralization  factor  ^  is  assumed  to  be  a 
constant.  This  is  a  reasonable  approximation  when  /  <  .5  .  For  /  >  .5, 
since  the  return  current  density  has  a  radial  profile  considerably  different 
from  Jbo>  the  r.h.s.  of  (28)  is  not  valid. 


Since  the  equilibrium  distribution  function  fQ  is  an  integral  of  motion, 
it  can  be  written  as  a  function  of  the  constants  of  motion.  We  choose 
the  separable  form 

f0  =  G j_(w)  F(v)  (30) 


where 


w  =  v 2  - 


2eZ.R  R 
b  oo 

Y  mvZ 
'  o  o 


(31) 


is  a  constant  of  motion  according  to  Eq.  (26)  and  the  velocity  spread  v 
is  assumed  constant  by  our  model  assumption  (8).  G  and  F  are  normalized 


such  that 


"o<r>  *  f<\  G1 

(32) 

1  =/dv  F 

(33) 

where 

nQ  is  the  equilibrium  particle  density. 

For  the  Maxwellian  distribution,  G  and  F  take  the  form 

G„  .  \  e 

HUj^ 

(34) 

F  _L_e-(vW) 

/2v  u 

(35) 

where 

the  mean-square  velocities  are  defined  by 

2  2 

U1  =  <vf> 

(36) 

2  2 

u  E  <v  > 

(37) 

The  normalization  constraints  (32)  and  (33)  are  appropriately  satisfied. 
The  transverse  distribution  GM  is  the  important  isothermal  case  described 
by  Bennett  (Ref.  8  ).  It  has  been  shown  by  Lee  et.  al.  (Refs.  4  and  9  ) 
that  the  Bennett  distribution  is  a  good  approximation  to  the  pinch- 
equilibrium  state. 


It  can  be  readily  verified  that  the  Bennett  profile 

Jt>°  =  °B  E  ^  a  -\W)Z 


and  the  magnetic  potential 


•  AB  -  -  '--'tr—  l"  11  +  r2'*2) 


(38) 


(39) 
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satisfy  Eqs.  (28),  (29),  (31),  arid  (34),  with 


I 


b 


eZ,  n 
b  o 


(40) 


Taking  the  virial  moment  of  (26),  the  mean-square  transverse  velocity 
satisfies  the  equilibrium  condition 


2  2 
<kB  r  > 


=  (1  -  tm) 


(41) 


with  the  Alfven  current  1^  defined  by  (1).  This  velocity-spread  is  due 
entirely  to  betatron  oscillations  --  energy-spread  effects  are  neglected. 


For  later  reference,  we  express  the  longitudinal  mean-square  velocity 

2  2 
spread  u  in  terms  of  the  energy  spread  t£  =  (y  -  y  )mc  .  From  Eqs.  (13) 

through  (15),  we  derive  the  relation 


2  2 

u  _  <v  > 


1 <  ( AE)S 


i  2 

(y0  -  1) 


Yomc 


(42) 

(43) 


Note  that  if  we  attribute  the  longitudinal  veloci ty-spread  to  betatron 
motion  only.--i.e.,  ignoring  energy-spread-'then 


jbetatron 


-4 


0  0 


rj  2  1. 

2y0  A 


(44) 


which  places  a  lower  bound  on  the  longitudinal  temperature. 


As  a  consequence  of  ignoring  the  v-dependence  in  k^(r),  Eq.  (27),  the  radial 
profile  of  the  equilibrium  current  is  independent  of  the  particular  subset 
of  beam  particles  with  the  same  v.  This  is  expressed  explicitly  by  intro- 

A 

ducing  the  partial  current-density  for  the  subset  of  particles 


with  longitudinal  velocity  v  =  v^  +  v  :  i.e., 

3b0(r.v)  =  r(v)  Jb0(r)  (45) 

Jbo(r)  =  /dv  Jb0(r-v)  (,|6) 


In  subsequent  sections,  we  shall  treat  each  subset  of  particles  separately 
in  order  to  account  for  the  longitudinal  velocity  effects  on  hose. 
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IV. 


HOSE  DISPERSION  EQUATIONS 


Hose  instability  is  treated  by  considering  small  amplitude  perturbations  of 
the  form 


6  J  i  =  J ,  ~  J ,  — 

b  b  bo 

,  JO 

Jbl  e 

(47a) 

6Jb(v)  =  Jb(v)  - 

^bolv)  =  ^bl(v)  e° 

(47b) 

6A  =  A  -  A  =  A. 

o  1 

e10 

(47c) 

6f  =  f  -  fo  =  fi 

ei0 

(47d) 

where  the  zero  subscript  refers  to  the  axi -symmetric  equilibrium  state 

A 

discussed  in  the  last  section.  J^(r,' ,z,v)dv  represents  the  partial  current- 
density  of  beam  particles  with  longitudinal  velocity  between  (v  +  v)  and 

A  U 

(v  +  v  +  dv).  The  perturbed  amplitudes  (d,.,  A.,  f . )  are  assumed  to  be 

1  in  • 

small.  Note  that  the  azimuthal  dependence  e  in  (47)  is  selected 
specifically  for  the  hose  mode. 

In  this  present  study,  we  consider  the  longitudinal  velocity-spread  effects 
in  the  contexts  of  two  models:  the  rigid-beam  and  field  model;  and  the  multi- 
disk  distributed-mass  model  of  Reference  7. 
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IV. A.  Rigid-Beam  and  Field  Model 

Following  the  rigid-beam  development  of  Reference  7  ,  we  represent  the 
perturbed,  partial  current-density  to  lowest  order  by 


;bl 


dJ,„ 

=  'Xv<^  W 


(48) 


An  equation  of  motion  for  the  transverse  displacement  amplitude  can  be 
derived  by  averaging  Eq.  (23)  over  all  particles  in  the  disk  (r,,z)  with 
v  (=  vq  t  v).  Employing  the  rigid-field  assumption, 

dA„ 

Aj(C)  =  -  0(U  Tr-  (49) 


we  find 


T5TVT7 

o  o  b 


/ 


rdr  ^boY 


(50) 


With  the  Bennett  profile  JbQ  =  Jg  ,  Eq.  (38),  (50)  simplifies  to 


with 


d  .X 


dz 


2“  =  kjB  '  XV^>Z)J 


(51) 


4 5(1  •  o  1 


eZ.R 
b  o 


m  d  6Tia^ 

o  o 


2  ^  '  U  !b 

"  3  a2  *A 

The  perturbed  form  of  Ampere's  law  (22)  is 


(52) 


(53) 
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d_/l  d_ 
dr  \r  dr 


Rrtov„ 
o  o 


'Jbl  =  'A  ^bl 


(54) 


Substituting  the  rigid-beam  and  field  representations  (48)  and  (49)  into 
(54),  multiplying  by  dAQ/dr,  and  integrating  over  the  radius,  we  find  the 
simple  relation 


where 


D(£)  ,  Ti  ^ 


<xU,0>  -/*  F(v)  \(C,z) 


(55) 


(56) 


For  the  Bennett  equilibrium,  the  dipole  magnetic-decay  length  T.  is  defined 

by 


3R  B  /* 00 

,  =  - j-  I  r  dr  o(r)  - 

1  2  a2  Jn  (1  + 


For, 


we  have 


2  a  */o 

o(r)  =  a 


(1  +  r2/a2)2 


0  (1  +  r2/a2)2 


(57) 


(58) 


T1  = 


RoV  6o 


(59) 


Equations  (51),  (55),  and  (56)  specify  completely  the  generalization  of  the 
rigid-beam  model  for  longitudinal  velocity-spread.  It  differs  from  the  zero 
temperature  model  (Ref.  7  )  by  the  full  z-derivative  in  (51)  and  the  v  - 
average  (56). 

If  initial  conditions  are  imposed  at  z  =  0  for  £  >  0  (where  the  pulse  head 
begins  at  £  =  0),  then  a  Laplace  transform  analysis  of  Eqs.  (51)  and  (55) 
results  in  growth  for  the  displacement  <X>  bounded  by 
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+<»  +  ic 


<X>  *  ~  f%m0)  e^o^> 

+io 


(60a) 


where  H(f?o)  depends  on  the  initial  perturbation  of  the  pulse.  The  growth- 
phase  exponent  g  is  expressed  by 


9  = 


(60b) 


where  u)(fi  )  is  extracted  from  the  single-mode  dispersion  equation  derived 
below.  Note  that  the  contour  of  integration  lies  iri  the  upper  (Ini  >  0) 
half-plane. 


In  a  single-mode  approach,  we  assume 


D,  Xv,  <x>  «  g-'Kkz-ujt) 


(61a) 


Transforming  to  the  independent  variables  (E,z),  Eq.  (20),  we  have 


D,  X  ,  <X> 


e-"(Vvo)z 


+  (w/vQ)£) 


(61b) 


where 


(l  =  u>  -  kv 
0  0 


(61c) 


Substituting  (61)  into  Eqs.  (51)  and  (55),  employing  the  z-derivative  (21c), 
we  arrive  at  the  rigid-beam  dispersion  equation 


-1“(no>»’  TT^-g-  /' 


dv 


EM 


ll-(no-kv)  3 


where  we  have  trarsformed  to  the  dimensionless  units 


(62) 


n 

TT~  ■*  fi0  •  IT"  ^  k 
o  SB  0  ksB 


(63a) 
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TlksB 


H  A 


(63b) 


v 


=  V  , 


We  now  consider  various  specific  forms  for  the  distribution  F(v). 


(63c) 


We  recover  the  cold  fluid  (or  zero  temperature)  results  of  Reference  7  by 
choosing 


Fc(v)  =  6(v)  (64) 

With  this  distribution,  (62)  becomes 

- - - j-  -  l  (65) 

(1  -<„)(!-  fl^l 

which  agrees  with  the  usual  form  of  the  rigid-beam  dispersion  relation  -- 
with  the  exception  of  the  current-neutral ization  factor  in  the  denominator 
which  is  usually  set  equal  to  zero. 


For  finite  temperature  hose,  the  Lorentzian  distribution 


Fl(v) 


A 

71 


2  2 
+  A 


(66) 


is  particularly  convenient  since  the  integral  in  (62)  can  be  performed 
analytically.  Any  function  K(v)  which  is  analytic  in  the  lower  half- 
plant  (Im  v  <  o)  with  peak  amplitude  less  than  |v|  satisfies  the  relation 


/ dv  Fl(v)  K(v)  »  K( v  =  -i  A) 


(67) 


Hence,  (62)  becomes 
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1 


(68) 


-M«Ja  = - - - =- 

U  -  yil  -  (%  +  i k A) ^ ] 

This  relation  is  analysed  in  Section  V  whore  its  consequences  on  hose 
growth  are  discussed. 

The  Lorentzian  tends  to  overestimate  the  effects  of  longitudinal  velocity 
spread  since  the  mean  of  is  divergent.  However,  in  the  limit  ■+  0, 
(68)  agrees  with  the  Maxwellian  result  [employing  the  distribution  (35)] 
to  lowest  order  in  ku  and  kA  if 

u  =  A  (69) 

Numerical  comparisons  suggest  that  the  real  part  of  the  r.h.s.  of  (68) 
agrees  reasonably  well  with  the  Maxwellian  weighted  results;  differences 
in  growth  rates  are  small  (<  40%). 


» 


) 


IV. B.  Distributed-Mass  Model 

The  distributed-mass  model  was  introduced  by  Lee  (Ref.  7  )  to  incorporate 
phase-maxing  --  associated  with  the  spread  in  betatron  frequency  -- 
into  the  rigid-beam  model.  Justification  for  the  model  relies  on  the  agree¬ 
ment  of  its  predictions  with  the  results  of  experiment  and  particle  simula¬ 
tions.  In  this  section,  we  generalize  it  to  account  for  longitudinal  velocity 
spread.  The  procedure  is  essentially  identical  to  the  rigid-beam  development 
in  Section  IV. A. 


Following  Reference  7  ,  we  associate  with  each  subset  of  particles  in 
(v,v+dv)  a  mass  distribution 


(70) 


where  n  is  a  continuous  variable  with  range  0<  y  <  1,  and  km  is  the  maximum 
betatron  frequency  for  a  particular  beam  profile,  i.e., 


s  <-■  ■  »> 


(71) 


From  (27)  and  (28),  we  find 


2  Jb0(°) 

vr  =  2  it  (1  -  {,  )  — 

m  '  Dnr  1^ 


which  for  the  Bennett  equilibrium  gives 

,2  _  2d  -  O  Jb 


(72) 


(73) 


A 

We  let  x  (C’z)  represent  the  displacement  of  a  disk  of  particles  at  (r,z) 

JJ  >  V  * 

with  mass  m^  and  velocity  v  =  vQ  +  y  The  center  of  current  is  now  given  by 


«x(s»z)» 


s  J“d\i  h  (y )  fc 


dv  F(v)  x  (£.z) 

M  >  v 
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where  h(y)  dy  denotes  the  fraction  of  disks  with  mass  m  and  is  normalized 

y 

to  one.  Following  the  arguments  of  Ref.  7  ,  we  choose 


h(y)  =  6 p  (1  -  y) 


which  is  zero  at  y  =  0,1  and  peaks  at  y  =  1/2. 


/s 

Substituting  (70)  into  (50),  we  get  a  new  equation  of  motion  for  X 


f‘  1D(«  '  xu.v(«'z)i 


and  Ampere's  relation  (55)  is  modified  to 


D(r)  +  T  -50.  =  >>  (77 ) 

These  two  equations  along  with  (74)  completely  specify  the  distributed-mass 
model  with  longitudinal  velocity  spread. 

The  single-mode  dispersion  equation  for  the  distributed-mass  model  follows 
substituting  (61)  in  Eqs.  (76)  and  (77): 

-MfiJA  *  fiv  Ml  -  u)  S  TwT  f dv - -  1  !  (78) 

0  Ja  (  11  <m)  L  v -(Vk'>)  ) 

where  the  dimensionless  units  now  correspond  to  setting  km^  =  1  rather 

than  k  D  in  (63).  Note  that  for  zero  velocity-spread  --  i.e.,  F(v)  =  F  (v) 

in  (64)  --  we  recover  the  distributed-mass  dispersion  result  of  Lee  with 

the  addition  of  the  current  neutralization  factor  /  . 

°nt 


A  tractable  solution  of  (73)  for  finite  velocity-spread  follows  by 
employing  the  Lorentzian  distribution  FL  of  (66).  Since  y  is  real  and  the 
singular  points  { kvg )  (from  setting  the  denominator  of  the  integrand  equal 
to  zero)  lie  in  the  upper  half-plane  (for  Im  >  0),  we  simply  replace 


kv  by  -ikA.  The  integration  over  p  can  now  be  performed  to  give 


-*“(", „>x  =  XT'“7mT  <4  +  G(fio  +  ikA» 


where 


G(S)  =  6S 
for  Ini  S  >  0  and  0  <  Re  S  <  1 


3  c2  .  c2n  c2v  Jl  -  S' 
2  -  S  +  S  (1  -  S  )  y.n^- — y- 


(79) 


(80a) 


The  branch  points  of  6(S)  lie  at  S  =  O,-1!  .  Since  the  9  -contour  of  inte¬ 
gration  in  (60)  lies  in  the  upper  half-plane,  we  choose  branch  lines  that 
extend  to  -i®  in  the  lower  half-plane.  On  the  primary  sheet  for  Ini  S  =  0, 
G(S)  is  analytically  continued  into  -1  <  Re  S  <  0  and  j Re  S|  >  1  as  follows 


+iir  +  £n 


Jtn 


1-S 

-S2 


-in  +  £n 

2 
9 

s 


1-S 

? 

1-S 

<2 


Jtn 


1-S 

2 


,  0  <  Re  S  <  1 
,  -1  <  Re  S  <  0 
,  | Re  S|  >  1 


(80b) 


For  Im  S  f  0,  the  prescription  for  determining  G(S)  is  now  completely  specified. 


In  the  limit  n  ->  0,  and  p  ->  0  (or  1),  the  Lorentzian  distribution  in  (78) 
reproduces  the  Maxwellian  result  [employing  the  distribution  in  (30)]  to 
lowest  order  if 

A  =  u  =  <v2>  (81) 


Generally,  however,  the  Lorentzian  result  (79)  by  equating  A  with  u  tends 
to  overestimate  the  effects  of  longitudinal  velocity-spread  --  but  by  less 
than  40%.  This  is  more  than  adequate  for  our  present  purposes. 


The  dispersion  equation  for  the  distributed-mass  model  with  finite  velocity- 
spread  is  now  completely  specified.  Its  implications  on  hose  growth  are 
explored  in  the  following  section. 


V. 


ANALYSIS  Of  DISPERSION  RELATIONS 


V.  A. 


Rigid- Beam 


Growth  in  the  beam  displacement  ■  as  a  function  of  ( f,z)  is  measured  by 
the  leading  amplitude  (60).  The  growth- phase  exponent  g  for  the  rigid-beam 
model  follows  from  (68);  with  k  g7.  >z  and  k  g<,  r,  we  have 


g(f0)  =  -iv  ‘ 

(ftn  +  ik/)2 
(I?JA  =  -  0 


-iw 


RBV  o 


I1  -  (%  +  ik/)  ] 


(82a) 

(82b) 


where  we  have  set  -  0  for  convenience.  For  the  case  k / .  =  0,  Figure  1 
gives  the  Ini  t,ug  =  0  contour  in  complex  ft  -space;  the  shaded  region  corresponds 
to  Im  >  0  i.e.,  where  hose  growth  occurs  --  and  the  broken-line  marks 

the  inversion  ^-contour  of  integration  for  Equation  (60).  The  ^-contour 
is  not  permitted  to  lie  entirely  in  the  1  owe r  half-plane  since  there  are 
intrinsic  singularities  at  ft  =  +1.  Since  1m  ft  and/or  Im  wDP  must  be  qreater 
than  zero  over  some  part  of  the  fto-contour,  unstable  hose-growth  results. 


For  kA  >  0,  the  real  axis  (Im  ftQ  =  0  in  Figure  1)  shifts  up  by  an  amount  k/> 


For  k/\>  .3536,  the  unstable  region  lies  entirely  in  the  lower  half-plane. 
Thus,  the  ^-contour  can  be  taken  along  the  real  axis,  and  hose  stability  is 
established.  However,  the  rins- velocity  spread  u  must  satisfy  the  condition 


ku  >  .36  (ksBv0) 


(83) 


which  is  unrealistic  in  practice. 


For  a  general  understanding  of  hose  growth,  we  pass  the  inversion  contour 
through  a  saddle  point.  The  saddle  point  ftQS  satisfies  the  condition 


1^0 /ft 


=  0  =  -iz  + 


2  (ft  +  ik  a) 

Oa 


os 


H  -  (ftos  +  ik  A)  ] 


(84) 


V 

* 
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For  zA  »  C  and  *.<  1,  we  find 


^os  ~  ^ 

g(s<0S) 


f.  ?  .1 

-  +1  i 

_ _ ^ 

Az  2 

A 2 

-  ik  A 


-  iz  +  i 


zr\J 


z  k  a  + 


(85) 

(86) 


The  damping  of  hose  growth  due  to  longitudinal  velocity  spread  is  now  explicit 
in  the  -zkA  term  of  (86). 


If  we  maximize  g  with  respect  to  z  with  i,  fixed,  we  find 


3max 


4xk  A 


^os  0 


(87a) 

(87b) 


This  maximum  occurs  at 


z 


o 


4Afc2A2 


(87c) 


Thus,  for  a  -*  0,  we  have  unlimited  growth  for  all  ^  as  a  function  of  z  -- 
i.e.,  an  absolute  instability  --  which  concurs  with  the  rigid-beam  result 
of  Reference  7. 


For  finite  velocity  spread,  the  hose  displacement  amplitude  peaks  at  each 

C-disk  with  maximum  exponent  g  and  then  begins  to  damp.  The  maximum  dis- 

max 

placement  occurs  at  larger  zQ  with  increasing  £  according  to  (87c).  Thus, 
the  hose  displacement  convects  back  into  the  pulse  (towards  larger  ^)  as  it 
grows  in  amplitude.  It  is  straight-forward  to  verify  that  the  convection 
velocity  is  given  by 

v  =  -  4Ak2A2  (88) 

C  u 

3 

If  we  allow  a  maximum  displacement  gain  of  10  ,  for  example,  the  maximum 
pulse  length  Lp3x  permitted  for  hose  stable  propagation  is 
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wm 

f 


Lpax  -  28  A  kA  (89) 

according  to  (87b).  If  we  consider  a  proton  beam  with 
Eq  -  nic2  =  50  MeV 

Ib  =  10  ka  (90a) 

a  =  .5  cm 

we  find  .1  <  A  <  1  and  k$^  ~  .051.  Condition’  (89)  then  becomes 

Lpax  <  (5.5  meters)  (90b) 

sBvo 

Thus,  longitudinal  velocity  spread  alone  does  not  ensure  useful  hose-stable 
pulses  without  imposing  severe  conditions  on  the  mis  velocity-spread. 

For  (5  >  0,  the  growth-phase  exponent  (82)  aquires  an  additional  growth 
component 


+  l 


X 


(91) 


which  increases  the  maximum  displacement  at  £.  Since  we  are  interested  in 
the  relative  importance  of  longitudinal  velocity-spread  on  the  hose  mode,  we 
ignore  this  plasma  return-current  effect  in  the  present  treatment.  Note 
again  that  when  the  denominator  in  (91)  becomes  less  than  .5,  the  present 
model  begins  to  loose  its  validity.  For  significant  current-neutralization, 
we  must  take  into  account  the  proper  radial-profile  of  conductivity  and 
its  dipole  perturbation. 
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V.B.  Distributed-Mass  Model 

The  growth-phase  exponent  for  the  distri buted-mass  model  is  given  by  (79) 
with  kmBz  ->z  and  kmBC  ->C: 

g(s20)  =  -ifyz  +  G(qq  +  i k a)  ^  (92) 

with  G  defined  by  (80)  and  ^  =  0.  For  kA  =  0,  the  inversion  9  -contour  for 
(60a)  is  plotted  in  Figure  2.  The  shaded  region  corresponds  to  Re  G  >  0  or 
unstable  hose  growth.  The  branch  cuts  are  designated  by  dark,  wavy  lines 
extending  to  -i«>.  Since  the  inversion  f20-contour  must  lie  above  the  branch 
points  fi  =  ±1,0,  it  can  not  be  deformed  to  lie  outside  the  region  where 
either  Im  >  0  or  Re  G  >  0.  Therefore,  hose  instability  cannot  be  avoided. 

For  finite  velocity-spread,  the  branch  points  and  the  Re  G  >  0  shaded  region 
shift  down  by  an  amount  - i k A.  This  effect  is  illustrated  in  Figure  3  where 
Re  G  is  plotted  as  a  function  Re  for  various  values  of  kA  (with  Im  =  0). 
With  kA  >  .183,  the  unstable  region  lies  entirely  in  the  lower  half-plane  and 
the  (^-contour  can  be  taken  along  the  real  axis.  Thus,  we  achieve  absolute 
hose  stability  with  an  rms  velocity-spread  u  satisfying 

ku  >  .183  (kmBvo)  (93) 

which  is  nearly  a  factor  of  two  less  than  rigid-beam  condition  (83)  with  no 
kg-mixing. 

We  perform  a  saddle-point  analysis  on  (92)  to  extract  additional  information 
about  the  hose  instability.  The  saddle  point  9qs  satisfies  the  condition 


The  values  of  UQS  have  been  computed  numerically  (assuming  Im  fios  =  0)  and 
are  plotted  in  Figure  4  (solid-line)  as  a  function  of  kA.  Values  of  the 
maximum  growth-exponent  Re  G(fios)  have  also  been  determined.  For  these 
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values  of  S2os  and  kA,  we  compute  the  saddle  point  condition 


C  3(  In:  G )  s 

D(Re  it)  , 

'  os 

which  represents  the  rate  at  which  the  peak  wave-packet  convents  hack  towards 
larger  £,.  Values  of  v  / A  are  also  plotted  in  Figure  4  (dashed- line)  in  ten.', 
of  kA. 

For  a  given  kA,  we  can  compute  the  upper  bound  on  hose  growth  for  a  particular 
£,-disk  by 

<*> 

When  the  maximum  is  achieved,  the  displacement  amplitude  commences  to  damp. 

By  expression  (95),  wc  see  that  the  hose  growth  saturates  at  larger  7.  for 

3 

increasing  f.  If  we  allow  a  maximum  gain  of  20  ,  the  pulse  length  is  limited 
by 


RTg^T 


in  units  of  k" 


For  a  proton  beam  with  parameters  (90a),  we  have  for  kA  -  .05,  for  example, 
kmB  ~  ,08^  and 

Lpax  <  2.0  meters  [for  ku  =  .05  (k^vj]  (98) 

which  is  ^7.3  times  longer  than  that  allowed  by  the  rigid-beam  model  (90b) 
with  no  kg-mixing.  For  kA  =  .1,  we  find 


Lp  £3.6  meters  [for  ku  *  .1  (kmBv0)] 
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To  appreciate  the  effects  of  longitudinal  velocity-spread  alone,  we  compare 
(98)  and  (99)  to  the  cold-fluid  result  (kA  -  0): 

Lp3x  <1.2  meters  [ku  =  0]  (100) 

which  is  a  factor  of  three  less  than  the  hose-stable  pulse  with  ku  -  .1  (k  gV  ) 

Note  that  (98)  -  (100)  are  approximate  upper-bounds  on  hose-stable  pulse 
length.  Inserting  the  proper  units  back  into  (97)  and  employing  the  Bennett 
results  (59)  and  (73),  we  find  Lpiax  is  proportional  to  beam  current.  A 
decrease  in  beam  current  due  to  longitudinal  pulse-spreading  would  reduce 
the  hose-stable  pulse  length. 
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VI. 


CONCLUSION 


In  Section  V,  the  distributed-mass  model  with  longitudinal  velocity-spread 
provides  approximate  conditions  for  hose-stable  propagation.  It  is  useful 
to  compare  these  results  with  realistic  collective-accelerator  projections 
and  with  the  requirements  for  longitudinal  mass  stability. 


At  the  exit-port  of  the  Auto-Resonant  Accelerator  (ARA),  the  maximum  spre< d 
in  energy  is  estimated  to  be 


1 [II 

E 


max 


6.3  R 


m 


V1  -  P0> 


1/2 


(161) 


where  m  is  the  electron  mass  and  m  is  the  ion  mass  (Ref.  ]0  ).  This  result 
is  due  predominantly  to  longitudinal  velocity-spread  with  ion  axial • mot k  n 
confined  by  the  cyclotron  wave  potential.  For  the  50  MeV  proton  hear.  (9Cu), 
the  ARA  would  produce 


(AE) 


max 

- «  .04 


(JO  2) 


The  associated  rms  velocity-spread,  according  to  relation  (42),  is  so.icwi ■■.■n 
in  the  range 


.1  <  —  <  .2  (103) 

%  vo  ^ 

Comparing  this  to  predictions  (98)  and  (99),  with  t  ^  .5  kg,  we  would  expect 
a  significant  increase  in  hose-stable  pulse-length  due  to  longitudinal 
velocity-spread  at  the  exit-port. 


After  the  beam  propagates  a  distance,  some  of  the  longitudinal  temperature 
is  converted  into  transverse  temperature  by  the  azimuthal  magnetic  self¬ 
field  of  the  beam.  A  lower-bound  on  longitudinal  temperature  may  be  found 
in  the  dynamics  of  the  longitudinal  mass  (or  bunching)  instability.  Sloan 


f 


et.  al .  (Ref.  6  )  have  show*  that  the  longitudinal  bunching  instability  in 
the  axi -symmetric  monopcle  mode  is  Landau  damped  when 


u 

v 


2 

2 
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.  2 
a.  t: 

'  o'o 


(104) 


In  fact,  this  condition  is  enforced  by  !he  instability  itself:  as  the 
bunching  mode  grows  in  amplitude,  the  ras  velocity- spread  incrc-osos  until 
(104)  is  satisfic'd.  For  the  50  Me V  proton  beam  example,  we  find 


.075 


(101 


which  is  reasonably  close  to  the  lower  bound  estimate  (103)  produced  at 
the  exit-port  of  an  ARA.  Thus,  the  bunching  instability  appears  to  main  tain 
the  longitudinal  temperature  required  for  improved  hose-stability. 

K  •  conclude  from  these  preliminary  results  that  the  lcw-Y  light -ion  beam  ir-jy 
have  an  advantage  over  a  comparable  electron  beam  insofar  as  the  resistive 
hose  instability  is  concerned.  The  degree  of  advantage  (if  any)  depends,  of 
course,  on  the  amount  of  velocity-spread  the  pulse  can  tolerate  --  spreading 
of  the  pulse  length  would  reduce  its  current  and  decrease  its  hose-stable 
length.  It  has  been  suggested  that  a  "self-trapped  mode"  may  exist  where 
the  beam  particles  are  confined  longitudinally  by  the  self-induced  electric 
fields  at  the  front  and  rear  of  the  pulse  (Ref.  5).  Verification  of  this 
self-trapped  state  will  require  a  numerical  study  of  the  dynamic  coupling  of 
longitudinal  and  transverse  degrees  of  freedom  [e.g.,  a  modified  L0GAP  code 
(Refs.  1  and  2)]  . 


A  more  precise  quantitative  estimate  of  this  hose  advantage  would  also 
require  a  comprehensive  numerical  study  of  hose  dynamics:  a  sel f-consi stent 
model  with  a  realistic  o-channel  and  with  the  transverse  and  longitudinal 
degrees  of  freedom  properly  coupled.  The  latter  requirement  would  permit 
some  of  the  longitudinal  temperature  to  be  converted  to  the  transverse  plane 
by  the  magnetic  self-field.  This  would  result  in  a  new  isothermal  equilibrium 
state  with  a  modified  beam-current  radial-profile. 


VII. 
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In  the  ren::-. in  dor  of  this  soot  ion,  \;o  do.  '  •  t!  ■ 
govern  in g  equations  fc>r  .1  o*v-  y  propagation  an d  C  i  sc**;-  :  %.!■<_- 
numerical  n.c  fhods  of  rolu.  ion.  Tin  ru  w  ;onn:i-:; !  1  :md 
ulgori  thm,  GO,  which  j  vo.-:  axinl-symme  tr  j  c  ,  :vol  • 

equations  in  real  time,  has  p  d  eor.si  d-ruvbl  o  l  es  l.  i  ug. 
Some  remits  are  p;x sen  u  d  .  The  long  i  t  ud  i  nn  ]  di  spf-usj  o-\ 
algorithm,  bXGDFP ,  is  si  ill  under  develop:::*  ut  r  til  j c  a  . 

to  revision.  The  radial  envelope  a.  I  gor  j  lb;n  is  only  slight  i 
modified  from  the  version  contained  in  II I  CAP. 

The  low- y  propagat ion  code ,  LOGAP ,  is  intend jd 
to  model  concurrently  the  root-mean-square  ( rms )  radial, 
•envelope  and  the  longitudinal  charge  distribution  of  a 
single-pulse  in  real  time.  To  make  the  problem  tractable , 
v/e  impose  several  sirnpl  '  '  ./.i  ng  assumptions.  First,  wo  ignor 
the  complications  of  dipole  hose  instability ,  and  higher 
order  instabilities,  by  assuming  axial  symmetry .  Second, 
wo  impose  a  modi 1 ied-Bonnett  radial  profile  on  the  beam 
current  density.  This  reduces  a  many-dimensional  problem 
to  one  involving  only  the  rms-rad ius  and  its  time  deri vativ 
Third,  we  treat  the  radial  and  longitudinal  degrees  of 
freedom  separately. 

We  consider  a  single-pulse  of  ions  moving  in 

A 

the  positive  z  direction.  At  time  t=0,  the  mean  velocity 

v  (t=0)  is  assumed  to  be  known.  At  some  time  t>0,  vo  focus 

z 

our  attention  on  some  thin  slice  of  beam  cut  normal  to  the: 

A  - 

direction  z  and  moving  with  some  mean  velocity  v  (z,t). 

_  z 

Clearly,  at  constant  t,  if  v  varies  with  z,  then  beam 

particles  will  move  into  neighboring  slices,  resulting 

in  a  shift  in  v  and  charge  density, 
z 


A  1 


Y.’c  )..  .!!  i  t  U1  D.C  i  '.I  O  f  pa  l"l  1  C  1  t  ' ; 

ic.iily  by  j  n  irodueinj'  t!i .?  coord i  na  v«*  t  \  •>  or. 


K  ~  p  cl: 
o 


(  A  .  I  ) 


*  ^  Po  rt“K 


\.ho}.v  c  is  the  speed  of  li  f,ht  a  r. d  f  <*.  is  in  :  ti  nl  i  '  f  ■ 

o 

mean  pulse  velocity  at  t-0.  The  coordinate  r  j the  y  •: 
tion  of  our  new  reference  frame  in  uni  ts  of  It  np.th .  i 


displ  acemont  var  j  abl  e  C  denote:';  a  parti  cult.  tr 


cl  ice  of  space  measured  from  the  or.i in  of  the  lr.  ’  i : 
the  negative  a  direction .  The  origin  1,-0  is  ch  »:;«.•:)  t*' 
correspond  to  Die  front  of  the  puslc  at  l/- 1  -0.  Note  t  hat 
(A.l)  is  a  Galilean  transformation  -  we  are  still  in  the 
laboratory  frame  of  ref evince. 

Me.  tra.n.s  f  ora  our  a,o\a.  rn  i  np  equations  o  f  m :>  1  J  < 
with  the  following;  replacements : 


(:■  ♦  o') 

\dC  dey 


(A.  2) 


D_  _ 

*  Di 


For  a  fluid  descx*iption ,  the  Eulerian  time  derivative 


transformation  is 


Boc  h  *  (Bo°  ■  u’  !r 


(A. 3) 


where  v  (£,£)  represents  the  mean  vz  in  the  slice  of  beam 

at  <c,o. 


The  i*adi;\3  beam  current  densri  fv .  ,  and  l  be  lv  r  run  on  t 

density  in  the  z -d i recti on ,  J  ,  ore  assured  10  bo  hnuv.r: 
as  iunct  ions  o  f  ( r ,  C ,  C  )  .  The  conduct i vi  ly  c  is  r.i  nor::  tod 
as  a  function  of  (  r ,  f  ,  C  )  frc::.  the  vai  u<*:;  of  }:  ,  ]'  J. 

j  ‘  '* .  L  . 

For  our  particular  anpl  j  cut  i  on  ,  v.v  use  the  f: i ...pi  i  •  led 
induced  a:ir  conductivity  model  F/.'U  ;.i)  (ho  r.  y  )  . 


Decause  of  the  scn.ei  live  dope.- rah  -nee  of  the 
fields  on  conductivity,  it  is  necessary  to  first  pv  evib 
formal  solutions  to  Eq.  A  .4.  Wo  integrate  along  t  ho 
" c )  1  a r a  c teri stic  1 i n o s  ’ 1 , 


?A(C) 


o-iy 

—if  -  <r'"  01  +  «A<S> 


<1+Po> 


*b>  +  cB(cb) 


(A. 5b) 


4'C(C)  =  (C- 


(A .5c) 


where 


D 9_  <^o>  9_ 

i>CA  3C  30  3? 


(A .Ca) 


D _ 

dcb 


(1+3o) 
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(A.Cb) 


D 9_  a 

D?c  9?  3£ 


(A  .6c) 


A  4 


After  .i  mposi  ng  -symmetry  and  t)jo  coord  i 

transformations  ( A .  1 )  wj  Lh  (A. 2),  Maxwel.l  ’  y  hjul.1 
become : 
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(l-no) 
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o  t. 


(E'  +  ];.  )-  -p 
1  V  P 


.  7  o 

E'  +’  7°—  ]■:• 
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(A.4r,) 
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(A. 4b) 
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where  we  have  introduced  the  unit  conversion 
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E' 


EV 


7, 


Z. 


7 


free^space  impedance 
376.74  (ohms) 


=  0 


(A. 4c) 
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General ly 


if 


D  _ 


F(0  +  c(0  F(C)  -  0(0 


(A.  7a) 


thor; 


F( C )  -  P(CA)  e 


~£  g(C')d V  f  c(C  )cJC 

e  +e  CA 


(A. 7b) 


Hence,  the  exponential  dependence  of  fine  fields  on  conduc¬ 
tivity  is  made  explicit. 


A  5 
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f 


V(G  refer  to  the  liner,  (A  .  t» )  nr  churn 
In  F.i  euro  A .  1  ,  we  illustrate  our  d.i  f  fenuie  ing  gri 
c  h  a  r  a  c  1  c  ‘list  i  e  lint  ■  s  o  f  1  n  t  o  g  v  at  ion.  I  f  w  o  kno  w 
of  the  fields  at  C^  for  all  (£,r),  then  the  value 
Cr,  "  Ac  are  found  by  in  t  ograting  upward  along 

I*  X. 


otcr  j  ret  i  cs . 
d  and  i  he 
the  valu<  s 
s  at 
t  hi.: 


characteristics.  The  relevant,  difference  equations  ornery  <" 
directly  from  (  A.  4 ) ,  using  the  formal  solution  (A.  7),  and 
are  expressed  in  terms  of  the  two  variables  (c,r).  The 
method  of  solution  is  now  essentially  identical  to  that 
of  CATHY  in  Rex.  3.  The  field  values  arc  updated 
explicitly  in  C  and  implicitly  in  the  radial  dimension. 


Our  treatment  of  the  (B)  and  (C)  character j sties 
has  two  phases  (see  Figure  A.l).  First,  at  the  start, 
when  the  field  profiles  are  evolving  rapidly,  the  C -stops 
are  small  and  the  characteristic,  linos  extend  all  the  way 
down  to  the  S^-line.  After  the  fields  have  come  to  some 


quasi-equilibrium  state,  we  may 

our  (B)  and  (C)  integrations  at 

In  other  words,  EB  ,  H,B  and  E,f 

r  4>  z 

the  old  and  new  values  at 


increase  Ac  by  commencing 
the  previous  C  slice, 
are  interpolated  between 


2 

The  conductivity  o  ( I — 1 , J )  in  Figure  A.l  is 
updated  according  to  BMCOND  (Ref.  3).  t'e  modify  its 
incorporation  into  GEM,  however,  by  integrating  along  the 
proper  character istic , Eq .  A . 5c ,  rather  than  along  the  line 
of  constant  C  (the  "frozen”  approximation). 


To  complete  our  description  of  GEM,  we  provide 
a  prescription  for  the  source  terms  J^^(r,^,C)  an  d  Jbr<r’^> 
The  radial  dependence  of  both  are  presumed  to  be  known. 


A6 


M 


For  J  ,  v.c  impose  the  :  nod J  fi  oh-Br  and 
radial  profile  j;.i  von  by 


Jb„0',E,O  « 

irR 


V  9  2 

(1  ”  "  /.«  c  ) 


1'  ( t,  ■  t  K  >  C  )  v  o  9 

B  ‘J  (l  +  r"74r 


(A . fa  ) 


whore  Dio  radial  cuiof  C  if;  r;c* !:  by  1  bo  reLaion 


F  =  r  R 
C  '"'R  B 


(A. 8b) 


The  Bonne  It  radius  satisfies,  by  dr  fin?,  lion  ,  the  relation 


W 


(A. So) 


with  the  rms-radi us  H  defined  bv 


R  = 


(A . 8d) 


Once  the  cutoff  coefficient  CD  is  chosen  (usually  according 

it 

to  the  size  and  resolution  of  the  radial  r.cch ) ,  the  Bennett 
coefficient  follows  directly  from  Fq.  A . Sc  and  A . SO 
The  normalization  coefficient  K  is  set  by  Die  condition 


;c 

J  Jba  AUr)  -  M  (C , C>  vz<5,?)  0-9) 

O 

and  is  dependent  only  on  C^.  The  beam  charge  per  unit  length 
of  pulse,  N  .  is  related  to  the  beam  current  I ^  by 


N(C,C) 


Ib(g,C) 

v,CC.C) 


(A. 10) 


A  8 


and  to  the  beam  charge  density  p,  by 

o 


C,p)  =  fc  Pl,(r,C,C)  d(',.r2)  (A.  11) 


Thus,  is  completely  specified  by  R,  N,  and  . 


The  radial  beam  current  J,  is  assumed  to  have 

br 


the  form  similar  to  Eq .  A.  8a 


Jbr(r»C.O  =  “a;  N  U,  C)  vr(r,C,r.) 
7rRB 


(l-r2/)>c2  )2 
(l+r2/K 2  )2 


(A  .  12a) 


where  we  have  assumed 


vr  (r , £  >  C )  = 


r 

R(f,C) 


(A. 32b) 


The  dot  in  R  denotes  the  total  Eulcrian  time  derivative 
of  the  rms-radius: 


R  = 


(A. 13) 


The  relation  (A. 12b)  is  a  direct  consequence  of  the 
self-similar  expansion  ansatz  employed  in  the  derivation 
of  the  rms-radial  envelope  equation  (Ref.  3). 


The  elec  t  ron. ague  t  :;.c  fields  in  GKM 
specified  as  functions  c  >  i  (q,c)  by  prov.bbi  ny  1! 


z 


Both  R 


ii 


re  so!u;  ’  (.1 


1 1:  <  radial 


.re  corny!' 
i  l ,  !<  a ! ,  d 

. '.-a- 1  Ci ’jr.- 


equation  of  tr.ol  .i  oij  derived  by  ].•.  e  and  Coopc  r  (Lei.  4  )  . 

With  si  i  ght  modi  fi cation  ,  our  adaptation  i  or  LOG.'.:'  i  r; 
identical  to  the  HI  GAP  forma!  i  an  discussed  in  Ref.  3. 

Vie  need  only  to  replace  the  total  time  derivatives  by  the 
transform  (A.  3).  The  resulting  equations  differ  from  tin..;, 
in  HIGAP  by  the  requi  ron.ont  that  we  integrate  along  the 
now  charao tori s tic 


tD(c)  “  —r~  <c-cu)  +  ?D(c,y  c-io 


and  that  the  betatron  frequency  ci-. 
replaced  by 


e 

yfil 


(b  vey) 


r 


(Eq.  3.5,  Ref.  3)  be 

(A. 15) 


where  0  H  —  ,  m  is  the  beam  particle  mass,  and  e  is  the 
electronic  charge.  The  envelope  equations  arc  solved 
simultaneously  using  4th  order  Runge-Kutta. 
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V.'e  must  now  ck-r  >  v<  a  set. 


*  .  ••  f  M 


for  the  nie; i n  chu  )*:;•••  per  unit  ]•  !.■;, uh  and  ' 


v  Adopt  j  ng  a  fir  ;d  ck -script  .in  n  ,  vc  ini  rob;.. 


distribution  function  f  who  re  f  ( r  .  ,  x  ,  'j! »  t  )  e  - ,  d  :  us 


defined  to  be  the  number  of  be:-.:::  pa  r  l  i  cl.es  wit!:  (r  ,  ,  p } 

23-  '  ' 


in  t)ie  phase  volume  dr  d  x.  dp.  The  sing!  e  pari  i  e?  e 


momentum  *p  is  related  to  the  single  part iej  velour  ly  v  by 


p  «  yiov 


(A-  16) 


Because  we  normalize  f  to  the  total  number  of  beam  part  ir  ?. <• 


in  the  single  pulse,  we  have 


o  (  r  |  ,  z  ,  t  ) 


3 

f  dp 


(A.  17:0 


and  for  the  charge  and  current  densities, 


V*i '  z*  t)  =  e  n 


(A. 17b) 


,  v. ,  t)  =  e  I  v  f  dp 


(A .17c) 


lVe  define  the  beam-charge  line  density  by 
N(z,t)  =  fpb  dr± 


(A. 18) 


Sm 


I  i'  Q  represents  some  vari  able  of  j  r.t  crest , 
it  is  convc  n:  ■ -n  i.  to  introduce  the  notation  <0>  and  Q 
where 


<Q> 


if* 


X  lip 


(A.  3.9:0 


and  the  radio]  aver a re 


Q 


l  f 


Q  Pb  drx 


(X, .  10b) 


In  the  development  to  follow,  wo  assume  that 


Q  =  <Q> 


(A . 20) 


The  governing  equations  of  motion  can  be 
derived  by  starting  with  the  J.orentz  invariant,  col  list  on  less 
Boltzmann  equation: 

*?-£  +  *  Vf  +  F  ♦?  f  =  0  (A.  21) 

St  yin  p 

where  the  Lorontz  force  F  is  just 

F  =  eF  +  e~v  x  B  (A. 22) 

Since  we  are  interested  in  various  momentum  moments  of  A. 21, 
we  let  Q (p-)  represent  some  function  of  "p.  After  multiplying 

3 

(A. 21)  by  Q  and  integrating  over  /dp  ,  we  have 

~r  n  <Q>  +  ?  •  n  <q"v>  -  n  <F  •  V  >  Q  =0  (A.  23) 

dt  p 


A12 


where  we*  have  used 


P 


F  -  0 


We  now  multiply  Kq.  A.  2d  bv  the  electros!  at  i  c  e!v  e 

2 

and  porronr.  the  radial  integration  fdr^  to  find 


feK  <;P  -  f 


<  Q  v  >  -  h<  T*  V  Q>  0 
z  P 


We  have  imposed  a >; i  a  1  - ; :y  mm.c t r y  and  the  condition 


Lim  r  p,  -  0 

r-Hv> 


/■A  <■ i  r  ■  \ 

(P  .  2e  ) 


Finally,  by  performing  the  coordinate  transformation 
(A.  2),  we  have 


3C 


N  <Q>  + 


k  Ni<^  - 


<  F* 


V  Q> 

P 


r>  r*  \ 

2u  ) 


'  To  obtain  a  closed  set  of  equations  from 

Eq.  A. 26,  we  proceed  as  follows.  First,  we  substitute 
three  relations  for  Q: 


Q  -  1,  Pz,  and  y 


(A .27) 


where  p  =  ymv  .  For  our  particular  application,  we 
z  z 

treat  v  as  a  functional  of  (y,v.)  and  Taylor  expand  the 

Z  J- 

<">  terms  about  the  mean  gamma  y,  where 


Y  -  Y  +  Sy 

and 

<Fr>  =  0 


(A. 28a) 
(A. 28b) 
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Finally,  v:c  assume  7 ■■to  radial  correlations  between  v 
and  the  electric  field  Yt?,  i.e., 


<  v  y,  : 

z  z 


-  <Vz> 


:r:z> 


(A.  33) 


This  is  necessary  to  corv.pl  ote  the  separation  between  the 
transverse  and  longitudinal  mode .  The  mean  transverse 
velocity  is  assumed  to  be 

• 

*S3t 

v  =  (A.  34) 

in  keeping  with  the  self-similar  ansatz  of  the  radial 
envelope  equation. 
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(A. 35c) 
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|  <  e  E 

z>  +  <®VV’j 

(A . 35 e  ) 
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(A  .351  ) 
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z 

G2  K 

(A  .36) 
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V 

z 

0 2  v  K 
z 

(A  .37) 

The  source  terms  (A.35e,  A.  35f)  are  generated  in  GEM  using 
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reference 

frame  and  the  displacement  vaiui 
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,  c,  ,  denotes 

a  particular  transverse  slice  of  space  measured  iron  tin. 
or  join  of  the  moving  frame,  in  the  negative  z  direction. 

The  $  c  frame  velocity  .is  initialized  to  equal  the1  mean 
v  -velocity  of  the  beam  particles  at  t--0.  Note  that  these 
coordinates  are  still  in  the  laboratory  frame  of  reference 
The  resulting  difference  equations  follow  directly  from 
Eqs.  (A . d )  of  Aop .  A  after  imposing  the  formal  solution 
(A. 7).  The  "characteristic  lines"  and  the  interpolation 
scheme  are  also  described  in  Appendix  A. 
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is  violated.  (The  root  finding  algorithm  permits  a  con- 
vergence  error  of  1  0_  .  )  Our  particular  choice  for  Ar 
and  AS  is  set  by  the  dimensional  requirements  of  our  pulse 
profile.  Note  that  we  restrict  the  range  of  k^Ar  and  k^AS 
above,  rather  than  go  from  0  to  2rr .  This  is  sufficient  since 
the  remaining  range  would  give  the  complex  conjugate  equation 
where  |  9  |  =  |  fi  |  . 
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